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We study separability properties in a 5-dimensional set of 
states of quantum systems composed of three subsystems of 
equal but arbitrary finite Hilbert space dimension. These are 
the states, which can be written as linear combinations of per- 
mutation operators, or, equivalently, commute with unitaries 
of the form U ® U ®U. We compute explicitly the follow- 
ing subsets and their extreme points: (1) triseparable states, 
which are convex combinations of triple tensor products, (2) 
biseparable states, which are separable for a twofold partition 
of the system, and (3) states with positive partial transpose 
with respect to such a partition. Tripartite entanglement is 
investigated in terms of the relative entropy of tripartite en- 
tanglement and of the trace norm. 
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I. INTRODUCTION 

One of the difficulties in the theory of entanglement 
is that state spaces are usually fairly high dimensional 
convex sets. Therefore, to explore in detail the potential 
of entangled states one often has to rely on lower dimen- 
sional "laboratories" . An example of this was the role 
played by a one-dimensional family of bipartite states , 
which has come to be known as "Werner states" . In this 
paper we present a similar laboratory, designed for the 
study of entanglement between three subsystems. The 
basic idea is rather similar to Q, and we believe this 
set shares many of the virtues with its bipartite coun- 
terpart. Firstly, the states have an explicit parameteri- 
zation as linear combinations of permutation operators. 
This is helpful for explicit computations. Secondly, there 
is a "twirl" operation which brings an arbitrary tripar- 
tite state to this special subset. This proved to be very 
helpful for the discussion of entanglement distillation of 
bipartite entanglement: the first useful distillation pro- 
cedures worked by starting with Werner states, applying 
a suitable distillation operation, and then the twirl pro- 
jection to come back to the simple and well understood 
subset, thus allowing iteration |^,||. Geometrically this 
means that the subset we investigate is both a section 
of the state space by a plane and the image of the state 
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space under a projection. The basic technique for get- 
ting such subsets is averaging over a symmetry group of 
the entire state space. Such an averaging projection pre- 
serves separability if it is an average only over local (fac- 
torizing) unitaries. Of course, special subgroups might 
turn out to be useful. For example, in a recent paper Q| 
a class of tripartite {n = 3) states was studied for dimen- 
sion d = 2, which is invariant under unitaries of the group 
of order 24 generated by <rf 3 , cr 3 ® 1 <g> 03 , 11 C*D 03 03 
and exp(i7ra- 3 /3) 03 . 

The third useful property of the states we study is that 
they can be defined for systems of arbitrary finite Hilbert 
space dimension d, leading to the same 5-dimensional 
convex set for every d. (This generalizes to an (n! — 1)- 
dimensional set for n-partite systems). Surprisingly, it 
turns out that the separability sets we investigate are 
also independent of dimension. 

We now describe the natural entanglement (or sepa- 
rability) properties we will chart for these special states. 
Our classification is similar to the one used in B, but 
differs in that we do not artificially make the classes dis- 
joint. 

Of course, we can split the system into just two sub- 
systems and apply the usual separability/entanglement 
distinctions. A split 1 1 23 then corresponds to the group- 
ing of the Hilbert space 'Hi®TL2®'H3 into H.i®(1i.2®'Hz). 
We call a density operator p on this Hilbert space 1 1 23 - 
separable, or just biseparable if the partition is clear from 
the context, if we can write 

p = £a qP (W 23) , (i) 

a 

(23) 

with X a > and density operators p a on H2 <S> H.3. 
We will denote the set of such p by B\. This set will be 
computed in Section ^. Furthermore, as it is a necessary 
condition for biseparability (cf. Peres we are going to 
look at those states p having a positive partial transpose 
with regard to such a split, denoted by p G V\. Recall 
that the partial transpose A 1— > A Tl of operators on Hi (8 
H2 is defined by 

(£A,®5 Q ) Tl =Y^Al^B a , (2) 

a a 

where A T on the right hand side is the ordinary transpo- 
sition of matrices with respect to a fixed basis. It is clear 
that B\ C V\ holds, but as we will show in Section |VT|by 
computing Vi, this inclusion is strict except for d = 2. 
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As a genuinely "tripartite" notion of separability, we 
consider states, called triseparable (or "three-way classi- 
cally correlated" ) , which can be decomposed as 

p = E A «^WW 3) , (3) 

a 

where A Q > 0, and the p$ are density operators on the 
respective Hilbert spaces. The set of such density oper- 
ators will be denoted by T . Of course, we may also con- 
sider states which are biseparable for all three partitions. 
It is known |(| that this does not imply triseparability, 
i.e. T 5 (Bi n^fl Bz). Further examples will be found 
below. 

Since in this paper we will only be interested in a five 
dimensional set W of symmetric states (see the next sec- 
tion), we will from now on use the symbols T, B\ and V\ 
only for the corresponding subsets of W. 

II. DEFINITION AND MAIN RESULTS 

A. U ® U ® {/-invariant states: W 

Throughout we consider states on a Hilbert space of 
the form T-t®Tt®Tt, where TL is a Hilbert space of finite 
dimension d. The group of permutations on 3 elements 
acts on this space by unitary operators V n , defined by 

For the six permutations 7r we use cycle notation, so that 
V(i2) is the permutation operator of the first two factors, 
and V( 12 3) is the cyclic permutation taking 1 to 2. We 
denote by u dU" the normalized Haar measure on the uni- 
tary group of TL, and define on the space of operators the 
operator 

Pp = I dU (U ®U ®U)p{U ®U ®U)*. (4) 

Clearly, P takes positive operators to positive operators 
(it is even completely positive), and tr(Pp) = tr(p), i.e., 
P maps density operators to density operators. We can 
now define the set of states, which form the object of our 
investigation. 

Lemma 1 For an operator p onTL®TL®TL the following 
conditions are equivalent: 

1. (U ®U (g) U)p = p(U ®U ®U) for all unitary op- 
erators U onTL. 

2. Pp = p. 

3. p = ^2 P-kV-k with coefficients p^ G C. 

The set of density operators satisfying these conditions 
will be denoted by W. 



The equivalence of (1) and (2) is straightforward from 
the invariance of the Haar measure. The implication 
(3)=>(1) is trivial, because the permutation operators 
clearly commute with operators of the form (U ®U ®U). 
The only non-trivial part is thus (1)=4>(3) which is, how- 
ever, a standard result ( Q chap. IV) from representation 
theory. Of course, all this works for any number of tensor 
factors. 

The above Lemma does not address the question how 
to recognize density matrices in terms of the six coeffi- 
cients p^. Hcrmiticity requires p^-i = ~p^, leaving effec- 
tively six real parameters. One more is fixed by normal- 
ization, so that W is embedded in a five dimensional real 
vector space. In terms of the parameters p^ positivity is 
not easy to see. In order to get a better criterion it is 
best to study the algebra of operators, which are linear 
combinations of the permutations. The product of such 
operators can readily be computed by using only the mul- 
tiplication law for permutations. The abstract algebra of 
formal linear combinations of group elements (known as 
the group algebra) can be decomposed in terms of the 
irreducible representations of the underlying group, sug- 
gesting a basis which is much more handy for deciding 
positivity. Again this step works for any number of fac- 
tors, but we carry it out only in the case n = 3: We 
introduce the following linear combinations of permuta- 
tion operators: 

R+=-(t + V {12) + V(23) + V (31) + V (123) + F(32i)) (5a) 
R-=-(l- V {12) - V m - V (31) + V (123) + F (3 2i)) (5b) 
Ro = 1(2- 1-V (ia3) -V(32i)) (5c) 
Ri = l(2V m -V {31) -V m ) (5d) 

R2 = 7% ( V{12) ~ V{31) ) (5e) 

#3 = -^ (V (1 23)-%n))- (5f) 

Then i? + , i?_ , Rq are orthogonal projections adding up to 
1, and commute with all permutations. This means that 
they correspond to the irreducible representations of the 
permutation group: R+ and i?_ correspond to the two 
one-dimensional representations (trivial and alternating 
representation respectively), and these operators are in- 
deed just the orthogonal projections onto the symmetric 
and anti-symmetric subspaces of Tt (8 H. <S> H in the usual 
sense. Their complement Rq corresponds to a two dimen- 
sional representation, which is hence isomorphic to the 
2 x 2-matrices. The operators R\, R 2 , R3 act as the Pauli 
matrices of this representation. In other words, the six 
hermitian operators Rq, R 2 , R3 are charac- 

terized by the commutation relations RiR± = R±Ri — 0, 
Rf — Rq, for i — 0,1, 2, 3, and R\R 2 — iRz with cyclic 
permutations. 



2 



Now every operator p in the linear span of the per- 
mutations can be decomposed into the orthogonal parts 
R+p, R-P, and Rop, and positivity of p is equivalent to 
the positivity of all three operators. This leads to the 
following Lemma: 

Lemma 2 For any operator p onTl®Ti®Ti, define the 
six parameters rk(p) = tr{pR k ), for k G {+, — , 0, 1, 2, 3}. 
Then r k (Pp) = r k (p) . Moreover, each p € W is uniquely 
characterized by the tuple (r+, r_, ro, r±, r^, T%) € K 6 , and 
smc/i a tuple belongs to a density matrix p £ W if and only 
if 



E 



CkRk 



(J) 



r+,r-,r > 0, 
and 



r+ + r_ + r 



1 



r? < fn. 



(6) 



Note that in this parameterization the set W does not 
depend on the dimension d with one exception: for d = 2 
the anti-symmetric projection i?_ is simply zero, so for 
qubits we get the additional constraint r_ = 0. If one 
considers a given density operator p as an operator p' in 
H! ®TL' ® TL' for a higher dimensional space TL 1 D H, by 
setting all "new" matrix elements equal to zero, we will 
have r k (p) = r k (p'). 

Taking ro = 1 — r + — r_ to be redundant, we get 
a simple representation of W as a convex set in 5 di- 
mensions. Unfortunately, 5 dimensional sets are still not 
very amenable to graphical representation. For visualiz- 
ing the sets we are going to describe analytically, we will 
therefore use suitable 2 and 3 dimensional representa- 
tions. Again, we have the possibility of using sections or 
projections of W, and we will emphasize sections which 
can also be understood as projections. 

The simplest example of this is to take the subset 
W p C W of states, which also commute with all per- 
mutations. The corresponding projection is simply aver- 
aging with respect to permutations. Clearly, W p consists 
of those operators in W, which are linear combinations of 
i? + ,i?_,i?o alone. Takingr + and r_ as coordinates we 
get the triangle in Figure |l|. Thus each point in this tri- 
angle represents a density operator in W p . On the other 
hand, it represents the set of states in W projecting to it 
on permutation averaging: this will be all states with the 
given values of r + and r_ in the 6-tuple, which therefore 
differ only in the values of ri,r2, and r$. Thus over ev- 
ery point of the triangle in Figure [l] we should imagine a 
Bloch sphere of radius ro. 

If more detail is required, we will also use three dimen- 
sional sections and/or projections of a similar nature. For 
example, if we average only over the permutations V(23) , 
we get the subset C W with r 2 = r 3 = (see the 

dotted tetrahedron in Figure 10). Averaging only over 
cyclic permutations, we get the subset W cyc C W with 
r% = r2 = (which gives the same tetrahedron as W^ 23 ' 
with r\ substituted by r3.). 

We note for later use that the expectation values r k 
arc not the coefficients in the sum 



fc=+,-,0,l,2,3 



These are related to the r k by the following dimension 
dependent transformation (which is obtained by observ- 
ing that trl = d 3 , tr(V (12 )) = d 2 , and tr(V( 123 )) = d). 



r. 




FIG. 1. Description of W in terms of the triangle W 
the corresponding Bloch sphere for each point in W p . 



and 



r+= ^(d 2 +3d+2) c + 
r_= ^(d 2 -3d +2) c_ 
r l = ^(d 2 -l)c l for i = 0, 1, 2, 3. 



(8a) 
(8b) 
(8c) 



B. Overview of Main Results 

An overview of the main results of this paper is given 
in Figure |[ To keep the picture as simple as possible, 
we have only depicted the set W p , i.e., the triangle in 
Figure [l| Naturally, this reduction does not allow the 
representation of our full results, i.e., the detailed struc- 
ture of the five-dimensional convex sets T, B\ , and V\ , 
which will be described in the corresponding sections. 
However, we found this diagram quite useful as a basic 
map for not losing our way in five dimensions, and hope 
it will similarly serve our readers. 

The shading in Figure || marks different separabil- 
ity properties, and the points labeled with capital let- 
ters arise by projecting pure states with special proper- 
ties with the twirl projection (0). Some of these points 
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(D,E,F) do not lie in the plane W p , i.e, they have non- 
zero coordinates {r-i,r%,r^). They are represented by 
white circles, in contrast to the black circles (A,B,C,G,H) 
representing permutation invariant states in the plane 
W p . 

The triseparable states correspond to the black triangle 
A (ABC). It is easy to see that any triseparable state pro- 
jected by permutation averaging to W p is again trisep- 
arable, i.e., the projection of T onto W p coincides with 
T n W p . The extreme points of this set are 



.4 
B 
C 



|123) — ► (1/6,1/6,0,0,0) 

|111) — ► (1,0,0,0,0) 

(| 1 1 1> - V3|112> + V3jl21) - 3|122))/4 

— > (1/4,0,0,0,0) , 



where the notation W — ► { r +i T ii r 2t ^3) indicates 
that the pure state |^)(^| is projected to this point 
by P from (Q). In other words, (\E f |i?fc\E f ) = r% for 
k = +,—,1,2,3. Note that all three vectors given are 
product vectors, the one for C being the product of three 
vectors in the "Mercedes star" configuration in the plane, 
at angle 120° from each other. 




74 ~ 1/2 ' ' 

FIG. 2. Subsets of W p with different separability proper- 
ties. Black: triseparable, dark grey: biseparable, light grey: 
images of biseparable states under permutation averaging. 
Special points explained in the text. 



A quantitative description of the genuinely tripartite 
entanglement of W is given in section [V in terms of the 
relative entropy and the trace norm. 

The biseparable set B\ is not permutation invariant, 



map B\ into itself, and we have to distinguish in our dia- 
gram between points (r + , r_) such that (r + , r_, 0, 0, 0) is 
biseparable (i.e., the intersection B\ D W p ), and points 
(r_|_, r_) such that for some suitable (rj., r2,ra) the quin- 
tuple (r_)-, r_, 7*1, r%, r^) represents a point in B\, (i.e., the 
projection of B\ onto W p ). In Figure || the intersection 
is the triangle A (GAB), drawn in a darker shade of grey 
than the triangle A(EFB), which is the projection of the 
biseparable subset B\ . Note that the shading reflects the 
inclusion relations, i.e., triseparable states are, in par- 
ticular, biseparable, and the section of the biseparable 
set is contained in its projection. Of course, the states 
in B\ n W p are also biseparable for the other two par- 
titions, since they are permutation invariant. Similarly, 
the projections of B2 and B3 onto W p are the same. 

Points of special interest for the biseparable set arise 
from the following vectors: 

-> (1/3,0,2/3,0,0) 
|121))/V2 — > (0,0,-1,0,0) 
|132»/v / 2 — > (0,1/3,-2/3,0,0) 
1 121) - v / 3|122))/v / 5 — ► (1/5,0,0,0,0) . 

Here the points B,D,E, and F are extreme points of f?i, 
and span a tetrahedron, which is equal to the subset 
B\ n W^ 23 ) of states invariant under the exchange 2^3. 
The point G lies on the line connecting E and D, and 
is the unique extreme point of B\ D W p which is not 
triseparable. In this sense it represents an extreme case 
demonstrating the inequality T 7^ (Bi n B2 H B3). 

The set V\ of states with positive partial transpose with 
respect to the partition 1 1 23 contains B\ strictly, but the 
difference cannot be seen in this diagram. In fact, we will 
show in Section [v| that even the 23-invariant subsets of 
V\ and B\ coincide, i.e., V\ H W^ 23 ) is spanned by the 
same four extreme points B,D,E, and F. 

As will be seen in Section VI there is a close connection 



D 


1122) 


E 


(|112) 


F 


(|123) 


G 


(|112) 



since the partition 1 1 23 clearly is not. As a consequence, 
the permutation average projecting W onto W p does not 



between the problems of finding V\ and finding states 
invariant under averaging over all unitaries of the form 
U ®U ®U . It turns out that the sets of triseparable and 
biseparable states commuting with such unitaries can be 
obtained via a simple linear transformation from their 
counterparts T DW and B\ n W computed in this paper. 
This mapping and a sketch of the results is given in the 
Appendix. 



III. TRISEPARABLE STATES: T 

If p is triseparable, hence has a decomposition of the 
form (||), we may also find a decomposition in which all 

factors pa are pure, simply by decomposing each of these 
density operators into pure ones. Applying to such a 
decomposition the projection P we find that p G T C W 
if and only if p is a convex combination of states of the 
form F* ( I \E r )(\E f | ) , where ^ = ip\ (g> -02 <8> "03 is a normalized 
product vector. Let us denote by 7^ ure C W the set 
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of such states. Our strategy for determining T will be 
to first get 7^ ure , and then to obtain T as its convex 
hull. The resulting characterization of T is formulated 
in Theorem ^. 

Given a product vector & = ipx ® tp 2 ® ^3> it is easy to 
compute the projected state P(|\E f )(\E f |): By Lemma|2]one 
just has to compute the expectations of the permutation 
operators. For example, 

(*|V(12)*) = (^1®^®^®^®^) = |(<0i|^ 2 )| 2 . 

In this way it is easily seen that the expectations of all 
permutations are {1, 01, a 2 , <i3, + io-5, «4 — ia>b\, where 
the 5 real parameters are given by 



fll = 


KV>2 |V3> I' 


(9a) 


a-2 = 


l(^3|^l)| 2 


(9b) 


«3 = 


l(^l|^2)| 2 


(9c) 


04 = 


S*5 «^l|^2><^a|^3><^3|^l» 


(9d) 


a 5 = 


^((VM^K^K^i))- 


(9e) 



Since a pure state in d dimensions (taken up to a factor) 
is given by 2d — 2 real parameters, these 5 quantities are 
a considerable reduction from the 6(d — 1) parameters 
determining the three vectors ipi. However, they are still 
not independent, due to the identity 



/(di, 02,03,04,05) := 04 + ■ 



aio 2 a 3 = 0. 



(10) 



Since we want to determine T pure exactly, we also have 
to find the exact range of these parameters, as the ipi 
vary over all unit vectors. This is done in the following 
Lemma. 

Lemma 3 A tuple (04,02,03,04,05) € R 5 arises via 
equations (^) from three unit vectors ipi, ip2,ip3 in o, d- 
dimensional Hilbert space (d > 3), if and only if equa- 
tion ( p7| ) is satisfied, < < 1 for i = 1, 2, 3, and 



1 — ai — a 2 — a>3 + 2a4 > 0. 



(11) 



If d = 2 the Lemma holds with last inequality replaced by 
equality. 

Proof: Necessity of equation (|To|), and < a, < 1 is clear. 
Inequality ( [Tl| ) is just the condition that the expectation 
of antisymmetric projection should be positive. Since 
this projection vanishes for d = 2, it is also clear that 
equality must hold in this case. 

Suppose now that ai,...,as satisfying these con- 
straints are given. We have to reconstruct ipi,ip 2 , and 
ipz satisfying equations (^|). These equations essentially 
determine the 3 x 3-matrix Mij = (ipi\ipj) of scalar prod- 
ucts. Of course, we already know the absolute values of 
its entries (note Ma = 1). The phases are irrelevant up 
to some extent: multiplying any row with a phase, and 
the corresponding column with its complex conjugate will 
not change the at after equation (^) , and amounts to mul- 
tiplying one of the ipi with a phase. Hence we may assume 



that the scalar products (ipi \1fj2) and (i/^IV's) are positive. 
The phase of the remaining scalar product (ipa\ipi) is then 
the same as the phase of 04 + 105, hence M is essentially 
uniquely determined by the <Zj. 

Now a matrix M is a matrix of scalar products if 
and only if it is positive definite: on the one hand, 
J2ijUi u jMij = || ~^2i Uiif>i\\ 2 > 0. On the other hand, we 
can construct a Hilbert space with such scalar products 
as the space of formal linear combinations of three vec- 
tors, with scalar products of basis vectors defined by M . 
Positive definitcness of M then ensures the positivity of 
the norm in this new Hilbert space. The dimension of this 
space is the rank of M (number of linearly independent 
rows/columns). So in the present case the dimension will 
be 3 (but any larger space will also contain appropriate 
vectors) or < 2, if M is a singular matrix. 

Positive definiteness of M is equivalent to the positiv- 
ity of all subdeterminants. The diagonal elements are 
1, hence positive anyway. Positivity of the three 2x2 
subdeterminants is equivalent to a.; < 1 for i = 1,2,3. 
Finally, the full determinant of M, expressed in terms of 
the Oi gives expression ([Tl|). It must be positive, and for 
d — 2 it must vanish, since M is singular. □ 



Lemma || describes the set 7^ urc of projected pure 
product states as a compact subset of the hypersurface 
in R 5 defined by equation (10). Computing the convex 
hull of this set in R 5 is the same as computing the con- 



vex hull of T^urc, because the expectations of permuta- 
tions or the operators Rk from (^) arc affine functions 
of the a,. Explicitly, the expectations = 
k = +, — , 0, 1, 2, 3, which we have used as our standard 
coordinates in W are 

r+ = - (1 + (ai + a 2 + <z 3 ) + 2o 4 ) 
6 

r_ = - (1 - (m + a 2 + a 3 ) + 2o 4 ) 
6 

ro = g(l - a 4 ) 
1 , 

n = - (2oi - a 2 - a 3 ) 
r 2 = -j=(a 3 - a 2 ) 

We begin by computing the projection of 7^, uro onto the 
(r+, r_)-plane, by determining the possible range of the 
combinations m = (04 + 02 + os)/3 and 04. By choosing 
phases for the scalar products we can make 04 vary in 
the range 1 04 1 < {aia2a 3 ) 1 / 2 — g 3 ^ 2 , where m and g are 
the arithmetic and the geometric mean of 04,02,03. As 
is well known, g < m, and equality holds if 04 = 02 = 03. 
Hence the projection of 7^ uro is contained between the 
parameterized lines 
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r + (to) = - (1 + 3to ± 2to 2/3 ) 
6 

r_ (m) = - (1 - 3m ± 2to 2/3 ) . 

Plotting these curves gives Figure [| It is clear that the 
shape is not convex, and its convex hull is the triangle 
ABC. 

T- 



0.2 n 

A 




FIG. 3. Section of the set T pule with W p and its convex 
hull. 



Proof: Suppose Xh is an extreme point of K. Then there 
must be a supporting hyperplane, i.e., a hyperplane H 
through Xh such that K lies entirely in one of the closed 
subspaces bounded by H . We claim that this implies 
that /, restricted to H, has to be either non-negative or 
non-positive in a neighborhood of Xh ■ 

Suppose to the contrary that there are points x+, x- £ 
HnU such that f(x+) > > We may then 

connect x+ and X- by a continuous curve lying entirely 
in U and also in one of the two open half spaces bounded 
by H. Since / is continuous, any such a curve must 
contain a point y with f(y) = 0, i.e., y £ (Nf f]U) C K. 
Since we can choose either side of H for the connection, 
we find points y £ K on both sides of H , hence H cannot 
be a supporting hyperplane. 

This argument shows, in the first instance, that the 
only possible supporting hyperplane at Xh is the tangent 
hyperplane (look at the Taylor approximation of / to 
first order). Applying the argument with the second or- 
der Taylor approximation, we find that hyperbolic points 
cannot have supporting hyperplanes, hence cannot be ex- 
tremal. □ 



A similar plot of the set 7p Urc including one more co- 
ordinate, r%, is given in Figure 0. 




1/6 1/4 1/2 r + 1 

FIG. 4. Plot of the same section as above making addi- 
tional use of the coordinate r^. 



Again it is clear that no point on the surface can be an 
extreme point of the convex hull of the surface, because 
the surface "curves the wrong way" . This is the intuition 
behind the following Lemma, by which we will show that 



also in the full five-dimensional case the interior of %_ 
contains no extreme points. 



pure 



Lemma 4 Let Nf = {x £ K" | f(x) = 0} be the zero 
surface of a function f £ C 2 (R n ,R), and K C E" a 
compact convex set. Let hi be an open ball around a point 
Xh £ Nf such that (U C\Nf) C K , and suppose that Xh is 
hyperbolic in the following sense: V f(xh) ^ ; and the 
tangent plane through Xh contains two lines such that the 
second derivative of f is strictly positive along one and 
strictly negative along the other. 
Then Xh is not an extreme point of K. 



To apply this Lemma to the function / from equa- 
tion (fL0|), we have to pick two appropriate tangent lines 
at any given point a = (ai, 02, 03, 04, 05) on the surface. 
We parameterize such lines as a + tb, t £ R so that 
f(a + tb) = /(a) + Alt 2 . Two choices with opposite 
sign of M are 



b = (0,0,0,05,-04), 
and b = (2ai, 2a2, 203, 3a4, 3as), 



M = {a\ + al) 
M = -3(a 2 + a 2 ), 



where we have used the equation f(a) = to evaluate 
the last expression. Hence every point of the surface Nf 
is hyperbolic. 

By Lemma ^|we therefore only have to consider bound- 
ary points of the surface, i.e., points for which at least 
one of the inequalities in Lemma |^ is equality. 

Let us begin with the equalities ai = 0, for at least 
one i £ {1,2,3}. Then we have 04 = 05 = by equation 



(TUT) and < aj + a/. < 1 (j ^ k) by equation (11). As 



we are looking for extremal points we are left with the 
cases aj = = representing the triorthogonal states 
B (i.e. point A=(|, i, 0, 0, 0) in the r^'s) or Oj +0^ = 1. 
All such points satisfy r_ = 0, hence they will be in our 
general discussion of cases with r_ = 0. The equalities 
Oj = 1 lead by ( |Ti"| ) to the inequality < 2a4 — (aj + o&) 
and therefore to 

1. , , , 

a-4 > -{aj + Ofc) > y/ajak = ^/alajak 



% > \ja\ = I a 4 1 > o 4 . 



From this we can see 05 = 0, 
Once again this implies r_ = 
only case to be checked. 



a j = a k aud 04 = aj = a^. 
so that this remains the 
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For r_ = 0, we can express the ai by r+, r%, r 2 , T3, and 
solve equation (|Io| ) for r3, obtaining a relation of the form 

r 3 =±h(r + ,n,r 2 ) , (12) 

where h is the square root of a third order polynomial. 
Equation (|l2] ) describes the surface of a convex set iff h 
is a concave function. This can be checked by verifying 
that the Hessian of h is everywhere negative semidefinite. 
Hence all points in T pnre with r_ = are extremal, and 
are characterized by equation This completes the 

determination of extreme points of T, summarized in the 
following Theorem. It also contains the dual description 
of T in terms of inequalities. 

Theorem 5 The subset T C W of triseparable states has 
the following extreme points, described here in terms of 
the expectations r^ = tx(pRk), k = +, — , 1, 2, 3: 

1. 3r| + (1 - 3r+) 2 = 

(r i + r + ) ■ (r i - \fir 2 - 2r+) • {r x + V3r 2 - 2r + ) 
and r_ = 0, 

2. The point A = (1/6, 1/6, 0, 0, 0). 

^4 siaie p£ W is triseparable if and only if it corresponds 
to the point A or the following inequalities are satisfied: 

(a) < r_ < \ 

(b) i(l - 2r_) < r+ < 1 - 5r_ 

(c) (3r 2 + [1 - 3r + - 3r_] 2 ) • (1 - 6r_) < 

(n + r+ - r_) • ((n - 2[r+ - r_]) 2 - 3r 2 ) . 

These inequalities are obtained by projecting the given 
point onto the hyperplane r_ = from point A, and to 
check whether the projected point satisfies the inequality 
l r 3| < h{ r +i r i> r 2) with h from equation ([l2|). To get 
an idea of the shape of T we compute the section with 
r+ = 0.27 and r_ = 0.1 (Figure |). 



r 3 




FIG. 5. Plotting the set T for the section r + = 0.27, 
r_ =0.1 gives a heart shaped surface with trigonal symmetry 
which is contained in the respective Bloch sphere. 



IV. RELATIVE ENTROPY OF TRIPARTITE 
ENTANGLEMENT 

Quantitative measures of bipartite entanglement and 
their properties are a very active area of research at the 
moment. In the tripartite case the difficulties in quantify- 
ing entanglement begin already with the pure states, for 
which no canonical form as simple as the Schmidt decom- 
position exists. One can, however, extend the standard 
definition of the relation "more entangled than" to tri- 
partite states. It is clear what local quantum operations 
should be in the multipartite case, and we can describe 
classical communication between many partners in much 
the same way as in the bipartite case. Once we fix the 
rules of classical communication (e.g., "each partner may 
broadcast her results to all the others" ) we will say that p 
is more entangled than a, whenever we can reach a from 
p by a sequence of local operations and classical commu- 
nication (LOCC), in which case we will write p >- a. 

A full characterization of this partial order relation is 
only known in the case of bipartite pure states (Nielsen's 
Theorem ||). Even in the mixed bipartite case there 
is no straightforward way of deciding whether one of two 
given density operators is more entangled than the other. 
Hence we cannot hope to give such a characterization in 
the tripartite case. Nevertheless, the entanglement or- 
dering is one of the features one would like to explore 
and to chart in W. There are two ways of approaching 
this: on the one hand, we may start from some state 
p € W, apply many LOCC operations to it, and see 
where we end up. We can always assume the operation 
to end up in W, because the twirl operation is itself a 
LOCC operation, which involves the random choice of U 
by any one of the partners, the broadcasting of U to the 
other two partners, and the unitary transformation by U 
at each of the sites. For an initial survey, we may even 
study the relation in the permutation invariant triangle 
W p even though the permutation of sites is definitely not 
a local operation. But if the inital state is permutation 
invariant, and T is any LOCC operation, involving cer- 
tain specified tasks for Alice, Bob and Charly, the three 
may just throw dice to decide who is to take which role. 
With this procedure they effectively get the permutation 
average of the output state of T. With such studies, we 
get sufficient conditions for p>- a. 

In order to get necessary conditions the only approach 
is to find functionals on the state space, which are mono- 
tone with respect to entanglement ordering. Luckily, one 
of the ideas for getting such monotones can be transferred 
from the bipartite case. Obviously, the triseparable sub- 
set is invariant under LOCC operations, so the distance 
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to T is an entanglement monotone, provided the distance 
functional has appropriate properties. One needs only 
one condition for a function A to define an appropriate 
"distance" A(p, a) between arbitrary states of the same 
tripartite system: 

A(Tp, Ta) < A(p, a) for any LOCC operation T . (13) 

Then for the functional 

E A (p)=M{A(p,a)\a€T} (14) 

we get the inequalities 

E A (Tp) < inf{A(Tp,cr)|c7 = Ta'; a' G T} (15) 
= M{A(Tp,Ta')\a' eT} (16) 
<M{A(p,a')\a' eT} = E A (p) . (17) 

Hence E A is indeed a decreasing functional with respect 
to the ordering y. Note that the only property of T 
needed to show this is that it is mapped into itself un- 
der LOCC operations. Any other set with that property 
(e.g., £>i or V\) will also lead to an entanglement mono- 
tone. 



Two natural choices for A satisfy requirement (13) 



and both of them satisfy it with respect to arbirtrary 
operations T (not just LOCC operations): firstly the 
trace norm distance: Ai (p, a) = \\p — <j\\\, and the rela- 
tive entropy As(p, a) — S{p, c), leading to entanglement 
monotones we denote by E\ and Es, respectively. In 
both cases, the actual computation of the distance for 
p, a G W is greatly simplified by the observation that we 
may consider both p and a as states (positive normalized 
linear functionals) on the algebra generated by the per- 
mutation operators, and that both the trace norm and 
the relative entropy are naturally defined for such func- 
tionals Moreover, because the twirl (||) is a condi- 
tional expectation the relative entropy of states in W is 
independent of the algebra over which it is computed (cf. 
Thm. 1.13, [0). Now the 6-dimensional algebra gener- 
ated by the permutations is independent of the dimension 
d, so that if we parameterize p and a by the expectations 
of Rk as before, we find that the entanglement mono- 
tones E A are independent of dimension. The expression 
for the relative entropy involves, apart from the abelian 
summands the logarithm of a 2 x 2-matrix, which can 
also be written explicitly in terms of the parameters rk 
for the two states involved. The variational problem flTi| ) 
can be then solved numerically for arbitrary states in W. 

The contour lines over W p of the resulting entangle- 
ment monotones are plotted in Figure ^ for E\ , and in 
Figure [?] for the relative entropy of tripartite entangle- 
ment Es- Note that the two neccessary conditions for 
p >~ a expressed in these diagrams complement each 
other. In order not to complicate these graphs we have 
not drawn the simplest sufficient condition for entangle- 
ment ordering: from any state p, any state lying on a 
straight line segment ending in T is less entangled than 
P- 



As a second section of interest we chose the plane 
r_ = = r\ — r2, which is relevant for qubit systems. 
Qualitatively, it gives the same picture of level lines wrap- 
ping around the tripartite set. 




0.2 0.4 0.6 0.8 

FIG. 6. Contour lines over W P for El 



r + 




FIG. 7. Contour lines over for Es. 



V. BISEPARABLE STATES: Bi 

In this section we are going to compute the set of bisep- 
arable states with respect to the partition 1|23. The tech- 
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nique is exactly the same as in the triseparable case: we 
first compute the set i3 puro of states of the form P*( | 1 ) 
with | X E')(5 , | biseparable, i.e., \& = ipi (g> ^2,3 ■ In a second 
step we get B\ as the convex hull of B puic . 

We are free to apply to our vector \t a U (8 U <S> U rota- 
tion without changing the projection. In this way we may 
choose ipi — |1). Now the rotated state *' is of the form 
*' = J2i j V'ijllu)- The expectations of permutations of 
such a vector, like 

j 

then depend linearly on the following real parameters: 

co = |V>n| 2 (18a) 
Cl | 2 (18b) 

c 2 = ^|^i| 2 (18c) 



C3 



c 4 + ic 5 



i>l 



(18d) 
(18e) 



From this we obtain the following r k : 



1 



r+ = - (1 + 5c + C! + c 2 + c 3 + 4c 4 ) 
6 



r*i 
f 2 = 

r"3 



(1 - c - Ci - c 2 - c 3 ) 
(1 - c - c 4 ) 
(-ex - c 2 + 2c 3 + 4c 4 ) 



3 

ci - c 2 

2C5 

V3' 



As in the tripartite case we need to determine the exact 
range of the parameters c$. Let us assume d > 2 for the 
moment. By the definitions of cq, c\ and c 2 we have 



co,c 4 ,c 2 > 0. 



(19) 



These parameters fix the weights of the blocks (i = 1, j = 
1), (i = l,j > 1), and (i > \,j = 1) in the normal- 
ization sum Y^i j—i l^ij] 2 = 1- c 4 + icg can be read 
as the scalar product of two (d — l)-dimcnsional vectors 
<Pi = ("012, ■ ■ ■ , i>id) and ip 2 = ($21, ■ ■ ■ , 4>di) with norm 



squares ||vi|| = C\ and ||(/? 2 | 
Schwarz inequality we have: 



c 2 . By the Cauchy- 



c 2 + 4 = KvMl 2 < IMI 2 M 2 = c lC2 , (20) 

and any value of c 4 + ic5 consistent with this can actually 
occur. 



We arrange the remaining ipij (i.j > 1) into a (ef— l) 2 - 
dimensional vector = (^ 22) . . . , ^ 2d; ^ 32j . . . , ifj dd ) with 
H^ll 2 = 1— c — ci— c 2 . On this (d—l) 2 -dimensional vector 
space, let U denote the operator swapping tpij and tpji. 
Then c 3 = (\I f |f/ 1 I f ) is the expectation of an hermitian 
operator with eigenvalues ±1. Hence 



M < ||*|| =1 -c - Cl 



C2 



(21) 



and all c 3 G R satisfying this inequality can occur. 

Together with the obvious modifications in the case 
d = 2, when there is only one index i > 1, we get the 
following Lemma: 



Lemma 6 A tuple (cq, ci, c 2 , c 3 , c 4 , C5) G R 6 arises via 
equations ( |7j| ) /rom a M77.it vector ^ in a d 2 -dimensional 
Hilbert space, if and only if equations jl^j), and (21) 
are satisfied, and, in the case d — 2, equality holds in (2L ) 
and (||7L). 



Let r denote the set of tuples (cq, Ci, c 2 , c 3 , c 4 , C5) sat- 
isfying these constraints. The depend linearly on the 
Ci, although the mapping is not one-to-one. Nevertheless 
any extreme point of B\ must be the image of an extreme 
point of the convex hull of T. 

Hence we can proceed by first determine the extreme 
points of T. Since the positive variables Co, |c 3 | and the 
sum (c\ + c 2 ) are only constrained by inequality j2l]), 
every point in T is a convex combination of tuples in 
which only one of these is equal to 1, and the other two 
vanish. This gives the extreme points 

1. c = 1 ^ r = (1,0,0,0,0) = B 

2. c 3 = +l^r=(i,0,|,0,0) = L> 

3. c 3 = -l^r= (0,|,-|,0,0) = F, 

and furthermore some points with (ci+c 2 ) = 1, Co = c 3 = 
0. Eliminating c 2 = 1 — c\ we can write inequality ( |20| ) 
as c\ + cl + (ci - 1/2) 2 < 1/4. This is a ball with extreme 
points parameterized by 



Co = 0, ci 



1 + cos(tf) 



C2 



cos(i?) 



sin(#) cos((p) sin(?9) sin(</j) 
c 3 = 0, c 4 = , c 5 = 



with ip, $ £ [0, 2tt]. By mapping this description of T to 
the r^-parameterization we come to the following Theo- 
rem: 



Theorem 7 The subset B\ C W of biseparable states 
with respect to the partition 1 1 23 has the following ex- 
treme points, described here in terms of the expectations 
r k = tv{pR k ), k = + ,-,1,2,3: 

1. The sphere given by j(3ri + l) 2 + 3r 2 + 3r| = 1 
with r_ = and r + = (n + l)/2 except for 



the point (|,0, ^,0,0), which is decomposable as 
(§,0 i 0,0) = |(S + J D) 
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2. The point F = (0, ±,-§,0,0) 

3. The point D = (±,0, |,0,0) 

4. The point B = (1,0,0,0,0). 

A state p G W is biseparable with respect to the partition 
1 1 23 if and only if it corresponds to the points F,B or D 
or the following inequalities are satisfied: 

(a) < r_ < i 

(b) -l < 1+r \z r 3r : 2r+ < 1 

(c) if -l < 1+r \z r i r : 2r+ < then 

3r%+3r 2 . + (l+2r 1 +r_-r + ) 2 < (2+r 1 -4r_-2r + ) 2 

(d ) if 0< 1+r \Z r s ~: 2r+ <1 then 

3r\ + 3r| + (1 - 3r_ - 3r + ) 2 < (r x + 2r_ - 2r + ) 2 . 

We omit here again the computation of these inequalities 
from the known extreme points. They can be obtained 
by projecting from the three points F, B and D onto the 
sphere of extremal points. 

The projection of the set Bi onto W p comes to be 
equal to the projection of the set of pure instates and 
was already shown in Figure || together with the section 
B\C\W P '. To compare B\ with T we plot again the section 
with r + = 0.27 and r_ = 0.1 (Figure g). 




embedded in the respective Bloch sphere together with T. 



To make the inclusion T C (B\ 062063) we mentioned 
in the introduction more evident we can now compute the 
sets £> 2 and B3 to build their intersection with B\. Due 
to the permutation symmetry of the three subsystems we 
can rotate B\ by ±^ in the ri-r2-plane instead. This 
leads to Figure |f 




FIG. 9. The intersection £>i ("162 063 is shown as mesh on a 
transparent surface allowing the set T to be seen. This plot is 
again computed for the section r+ — 0.27 and r_ = 0.1. The 
thick lines indicate the intersection of two of the biseparable 
sets. 



VI. POSITIVE PARTIAL TRANSPOSES: Vi 

One of the interesting aspects in the theory of bipartite 
entanglement to emerge in recent years was the consider- 
ation of the partial transpose of the density matrix, and 
in particular the positivity of the partial transpose. First, 
this positivity served as a necessary condition for separa- 
bility, which is even sufficient in 2(g) 2 and 2(g) 3 dimensions 
(the Peres criterion ||). Moreover, it is a necessary con- 
dition for undistillability, and here it comes much closer 
to sufficiency even in general situations. Both aspects 
play a role in the analysis of tripartite states. We will 
therefore describe in this section the subset V\ C W of 
states with positive 1-transpose. 

Since the dimensions for this bipartite system are 
d® d 2 , positive partial transpose does not automatically 
imply biseparability, i.e., the inclusion B\ C V\ may be 
strict. However, since we are considering a special class of 
states it is also possible that in this class equality holds. 
This does happen, for example, for the bipartite Werner 
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states . In the tripartite case we will see that B\ = V\ 
for d = 2, but not for higher dimensions, although the 
two sets come to be remarkably close (see Figure 13). 
However, the exact description of V\ is also important 
for distillation questions. 

The partial transpose A i— ► A Tl of operators on Hi <8> 
Ti.2 was defined in Equation (|^). In a tripartite system 
we take this operation to refer to the first of the three 
tensor factors, and write p S V\ if p Tl > 0. 



A. The algebra of partial transposes 

When p is a linear combination of permutation opera- 
tors as in Lemma [l], the partial transpose 



is likewise a linear combination of the six operators Vj 1 , 
and we have to decide for which coefficients /i w such an 
operator is positive. Since partial transposition is not a 
homomorphism, it would appear that the linear combi- 
nations of the Vj 1 can be a fairly arbitrary space of op- 
erators, and deciding positivity could be quite difficult. 
However, it turns out that these linear combinations do 
form an algebra, so after the introduction of the right 
basis deciding positivity is just as easy as determining 
the state space in Lemma |^. 

The abstract reason for this "happy coincidence" is 
that the operators Vj 1 span the set of fixed points of an 
averaging operation in much the same way as the per- 
mutations span the set of fixed points of P. The corre- 
sponding averaging operator is 



(22) 



p P = du (u®u®u)p(u®u®uy 



Its range consists of all operators commuting with all 
unitaries of the form U <S> U <£> U, hence is an algebra. 
The following Lemma describes the relation between P 
and P: 

Lemma 8 Let A be any hermitian operator, then 

1. PA = i« PA Tl = A Tl 

2. (PA) Tl =P(A T >) . 

Proof: For any hermitian operator A one has: 

PA = A & [U ® U <g> U, A] _ = 
<^ [U ®U ®U,A Tl }_ =0 
PA Tl = A Tl . 

Furthermore we can compute directly: 



PA 



''' - 1 dU (U®U® U)A Tl {U®U® uy 

dU {(U ®U ®U)A(U ®u ®uy) Tl 
= {PAf\ 



□ 



For deciding positivity of partial transposes we need 
a concrete form of the algebra spanned by the partial 
transposes of the permutation operators. For example, 
we get 



v (i2) = £(l*i*»ir = E \n k )(nk\ 



ijk 



ijk 



t, 



where $ = J^i I**) is a maximally entangled vector of 
norm d. The partial transposes of the other permutations 
are computed similarly. We can express all of them in 
terms of the first two: 



X 



V Tl 

V (12) 



and 



V 



V Tl 

"(23) 



V, 



(23) 



(23) 



as 



a Ti = i, v^ = vxv, v^ 23) = xv, v£ 21) = vx 



(321 



Then these operators satisfy the relations X* — X, and 
V* = V, and 



X v 



dX , 



V 2 



i 



XV X = X. 



(24) 



Due to these relations the set of linear combinations of 
the six operators {1, X, VXV, V, XV, VX} is closed un- 
der adjoints and products. Positivity of such linear com- 
binations, and hence the positivity of all partial trans- 
poses of operators in W can therefore be decided by 
studying the abstract algebra generated by two hermitian 
elements X and V satisfying (pij). As a six dimensional 
non-commutative C*-algebra it is isomorphic to the al- 
gebra generated by the permutations, i.e., a sum of two 
one dimensional and a two dimensional matrix algebra. 
But of course, the partial transpose operation mapping 
one into the other is not a homomorphism. 

From these considerations it is clear that all we have 
to do now is to find a basis of the algebra generated 
by X and V analogous to the basis (|g). This sort of 
computation can be quite painful, so we recommend the 
use of a symbolic algebra package. The result is 



s. 



a + v_r 2\ ^ ii + r 

2 



d+U 2 
2A \ 1 - V 



d-1 



S = -j^—r (d(X + VXV) - (XV + VX)) 



(25a) 
(25b) 
(25c) 
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Si = -^—r (d{XV + VX) -{X + VXvfj 



s 2 = 
s 3 = 



1 



i 



X - vxv 
XV -VX 



(25d) 
(25e) 

(25f) 



These operators satisfy exactly the same relations as the 
Rk from (|5|) and we will denote the corresponding expec- 
tation values by s&(p) := tr(pSk)- The two projections S± 
correspond to the two one-dimensional representations of 
the algebra, i.e., to the two realizations of the relations 
by c-numbers, namely X = 0, V = 1 and X = 0, V = — 1. 



can see in Figure 10. This is due to the fact, that the 
extremal points Pi and Qi (i = 1,2,3,4) lie on just two 
straight lines, namely P1Q4P4Q1 and Q 2 P 2 Q 3 Pj,- The 
intersection of the two tetrahedra is hence again a tetra- 
hedron, spanned by the extremal point s Po , Pi, Q3 and 
Qi (called E, B, F, and D in Sections II B and[V|), and 
is thus dimension independent. But it is easily verified 
from Theorem |?] that these four points are precisely the 
extreme points of the V(23) -invariant part of B\. Since 
B\ C Pi, we have shown the following: 

Lemma 9 A V(2 3 ) -invariant W -state has a positive par- 
tial transpose if and only if it is biseparable. 



B. The V(23) -invariant case 

The simplest case is the V^3) -invariant subset of 
W as it is a three dimensional object. In fact the 
V(23)-invariance implies the conditions tr(j(?V(23)) = lj 
tr(pV(i2)) = tr(pV (31) ) and ti(pV {123) ) = tr(pV" (32 i)). 
Therefore we have r 2 — r 3 — 0. In the same way we 
obtain for a V(2 3 ) -invariant state p S W Tl the condi- 
tions S2 = and S3 = 0. Positivity of a V(2 3 ) -invariant 
state in W requires now r + > 0, r_ > and |ri| < ro = 
1 — r + — r_ (cf. (^)) giving raise to a tetrahedron bounded 
by the hyperplanes 

(/ii) r+ = 0, (h 2 ) r_ - 0, (h 3 ) n = 1 - r+ - r_ 
(h 4 ) ri = r + + r_ - 1 

and having the extreme points Pi = (0,0,1), P 2 = 
(0,0,-1), P 3 = (0,1,0), and P A = (1,0,0). The same 
computation can be done on the partially transposed side 
leading to the tetrahedron confined by the hyperplanes 

(h[) S+ = 0, (h' 2 ) S_ = 0, (h' 3 ) 8i = l-8+- s_ 

(K) s i = s + + s - — i- 

Using Lemma ^|we can express the by the r/j of the cor- 
responding W-state. Multiplying by positive constants 
one gets an easier description of these hyperplanes: 

(h[) 2(1 + n - r_ - 2r+) + d(l + r x - r_ + r+) = 
(/i 2 ) 2(-l + ri + 2r_ + r+) + d(l - r x + r_ - r+) = 
(/l 3 ) 1 - r x - 5r_ - r+ = 
{h' A ) 1 + n - r_ - 5r+ = 0. 

Its four extremal points are now Qi = (^p, 0, ^jp), 
g 2 = (0, 2=S, -i±^), Q 3 = (0, i, - §) and Q 4 = (|,0, |). 
Of course, these point have no reason to correspond to 
positive states, and indeed only Q 3 and Q4 lie inside the 
state space, where Q\ and Qi are outside the state space 
for all d. 

As we are looking for those V(2 3 ) -invariant W-states 
that have positive partial transpose, i.e. that lie in Va, 
we have now to look at the intersection of these two tetra- 
hedra. The resulting object is again a tetrahedron as one 




FIG. 10. The two positivity tetrahedra bounded by the hi 
(dotted) and the h'i (dashed) and the intersection tetrahedron 
(solid lines) for d — 3. 



As we will see in the next subsection, the assumption 
of V(2 3 )-invariance is essential, i.e., the conclusion does 
not hold for general W-states. 

In order to see how V(2 3 )-invariance helps, we conclude 
this subsection with a direct proof of the above Lemma 
for d = 2. If p is a V( 23 ) -invariant W-state, then we can 
decompose it into the following sum 

P = 4(1 + %5))P(1 + V(23))+ \ (1 - V( 23 ))p(l - V(23)) 

=:p++p-. 
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It is now clear that p has a positive partial transpose 



According to this Lemma the set V\ can be visual- 



iff both p + and p~ each have a positive partial trans- ized as follows: firstly, one has to fix a point (r+,r_) 

in the permutation invariant triangle (Figure g). The 
possible choices of {r\, r<z, rs) can then be seen from Fig- 
ure |^. Apart from the heart shaped tripartite set in the 
center this Figure contains three quadratic surfaces: the 
Bloch sphere, and the two surfaces bounding B\. Com- 
paring condition (d) of Theorem [7] and the expression 
for Ri given in the above Lemma, we find that both 
constraints are given by the same hyperboloid, the one 



pose. p + denotes the V(23) -symmetric part of p, p the 
antisymmetric part. Thus we know that p + is a 2 x 3 
density operator and p~ a 2 X 1. For these systems the 
Peres criterion holds strictly i.e. states have a posi- 
tive partial transpose iff they are separable or in our case 
biseparable over the 1 1 23 split. Biseparability of p + and 
p~ is equivalent to the biseparability of p, which proves 
the lemma. □ 



C. The general case 

The positivity conditions for arbitrary linear combina- 
tions of the operators Sk give the following result: 

Lemma 10 Let p G W be a density operator with expec- 
tations r/j = tv{pRk), k = +,—, 1, 2, 3. Then the partial 
transpose of p with respect to the first tensor factor is 
positive, i.e. p £ V\, if and only if 



< r_ 

< r\ — r + — r_ + 1 
< 1 — ri — 5r_ — r + 
< -1 - n + r_ + 5r+ 

r| < R\ 

rl < R 2 



(26a) 
(26b) 
(26c) 
(26d) 
(26e) 
(26f) 



where 



i?i := (1 - n - 5r_ - r+)(-l - n + r_ + 5r + )/3 
i?2 := (1 — ri — r_ — r + )(l + r\ — r_ — r + ). 



Proof: Recall that averaging with respect to V( 2 3) projects 
V\ to the section of V\ with r 2 = ?'3. Therefore, the in- 
equalities describing the tetrahedron discussed in the last 
subsection are optimal. These are the first four inequal- 
ities. We therefore only have to describe the admissible 
set of (r2 , r3 ) , given (r + , r_ , ri ) . There are two conditions 
to consider, one from the positivity of p, and one from the 
positivity of p Tl . As shown in the first subsection, both 
these requirements have a very similar form, namely the 
positivity of an clement in an abstract algebra with two 
one-dimensional summands and one summand isomor- 
phic to the 2 x 2-matrices. Now in both cases (r + , r_, r\) 
are readily seen to fix the weights of the one-dimensional 
parts, as well as the trace and the expectation of the first 
Pauli matrix for the 2 x 2-part. This leaves a condition of 
the form r\ + r\ < R in both cases. The two conditions 
are given in the Lemma, where R2 = (1 — r + — r_) 2 — r% 
expresses the requirement p > 0. The condition ( |26q ) 
is obtained from p Tl > by expressing p Tl in the basis 
Sk, and applying the same criterion to the expectations 
s k - □ 



8L Hence in 
ing this hy- 



wrapped around the tripartite set in Figure 
that Figure we can readily find Vi by extenc 
perboloid all the way to the Bloch sphere, and taking the 
intersection. This is shown in Figure [ll], in the section 
r 3 =0. 



r 2 




FIG. 11. Plot of the Bloch sphere, 
r+ = 0.27, r_ = 0.1 and r 3 = 0. 



T, Bi and Vi for 



Figure |ll] shows the generic situation with r_ 7^ 0. 
When r_ = 0, in particular for systems of three qubits, 
the boundary ellipsoid of B\, described by condition (c) 
of Theorem [?], coalesces with the Bloch sphere. This leads 
to another instance where the Peres-Horodecki criterion 
for separability holds: 

Corollary 11 The intersections of B\ and V\ with the 
plane r_ — coincide. In particular, for 3-qubit W- 
states, biseparability is equivalent to the positivity of the 
partial transpose. 

We conclude this section by the explicit determina- 
tion of the extreme points of V\. 
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From Figure 

might appear that all points on the quadratic surfaces 
bounding V\ might be extremal. But this is misleading, 
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because we also have to take into account the possibil- 
ity of decompositions with different values of (r+ ,?*-). 
In fact, for the inequalities arising from p > it is 
evident that generically such decompositions are possi- 
ble: given any (r + , r_, r\, r%, rs), which lies on the Bloch 
sphere in Figure [Ly, we can just change the weights of the 
three blocks in the block decomposition of p according 
to (A+r+, A_r_, Aori, Aof2i Ao?"3)), as long as the \ a are 
positive, and the normalization A + r + +A_r_+A (l—r + — 
r_) = 1 is respected. This leaves a two dimensional afhne 
manifold through (r + , r_, r±, r 2 , r 3 ). Hence, unless other 
conditions constraining V\ prevent the indicated decom- 
positions no such point will be extremal. Of course, the 
second constraint (26e) has the same structure, because 
the algebra of partial transposes is isomorphic to the al- 
gebra generated by the states. Hence in Figure [n] only 
the points in the intersection of the hyperboloid and the 
Bloch sphere remain as candidates for extreme points. 
This is analogous to the extreme points of B\ , which also 
consist of the intersection of two quadratic surfaces in 
Figure 



11 



For V\ we get 



Theorem 12 The subset V\ C W of W -states with pos- 
itive 1-transpose has the following extreme points, de- 
scribed here in terms of the expectations r& = tv{pR^), 
k = +,-,1,2,3: 



1. The points P 2 , Q3, P4, and > 
V(23) -invariant part ofV\. 



}i, which also span the 



2. the remaining extreme points of B\, which form a 
sphere in the r_ = plane (cf. Theorem [?Jj. 

3. The points for which (r-j_, r_ , r\, 0, 0) lie in the in- 
terior of the V(23) -inv ariant tet rahedron, and for 
which inequalities (26e) and (26f) are both satisfied 
with equality. 

Proof: Let us first discuss the periphery of the tetrahe- 
dron. Every face of the tetrahedron corresponds to a face 
of "Pi, namely the face of points projecting to it upon 
V(23) -averaging. In Lemma [lO] this corresponds to the 
subsets for which one of the linear inequalities (26 a) to 
(p6|d) is equality. We will show first that each of these 
faces is actually contained in B\. Indeed, when (p6]b), 
(p6|c) or (p6|d) are equalities, one of the factors in R\ or 
R2 vanishes, forcing r 2 = r$ = 0, reducing our claim to 
Lemma|. When (p6|a) is equality, i.e., r_ = 0, the claim 
is contained in Corollary pT[ 

Now a point of V\ contained in one of these faces can 
only have decompositions in the same face, hence in B\, 
hence for such a point extremality in V\ and extremality 
in B\ are equivalent. 

It remains to show item 3 of the Theorem, i.e., to char- 
acterize the extreme points of V%, whose V/23) -averages 
fall in the interior of the tetrahedron. From the argu- 
ments preceding the Theorem it is clear that points for 
which only one of the inequalities (26c) and (261) are 
equalities cannot be extremal, since the surfaces defined 



by these equations contain straight lines. Therefore, the 
condition stated in the Theorem is necessary for a point 
to be extremal. It remains to show that none of the 
points with R\ — i? 2 can be decomposed in a proper 
convex combination. 

Let us denote by M\ (resp. M 2 ) the set of those points 
in the interior of the tetrahedron such that i?i < i? 2 
(resp. R 2 < Ri). The intersection M* = Mi n M 2 of 
these sets is described by the condition Ri = i? 2 , or 
explicitly 



r 1 + 3r_ + rir_ 



2rt 



3r, 



rir-| 



8r_ra 



2ri = l. 



(27) 



This is a one-sheet hyperboloid, generated by two sets 
of straight lines shown in Figure YL2. Consider a line 
segment 



(f + ,f-,h) + u(t+,t-,ti) 



(28) 



through one of the points p = (f + , r_ , f 1) € M* of the 
hyperboloid. Consider the radius functions y/Ri, evalu- 
ated as a function of the parameter u. If such a func- 
tion is affine (has vanishing second derivative) we can set 
(ri(w), 7"2(u)) = (cos a, sin a)\/R~i with arbitrary a, to 
get a straight line in the corresponding hypersurface in 5 
dimensions. We then call (t+,t-,ti) an affine direction 
for Ri. Along other directions, Ri is strictly concave, so 
no decomposition along the segment ( p8| ) is possible. For 
both radius functions, the set of affine directions is a two- 
dimensional plane, and thus best described by its normal 
vector. That is, t = (i_|_,t_,ii) is an affine direction for 
Ri if t ■ Ai = 0, where 



A l = 



2 - 3fi - 12f_ 
-2 - 3f 1 + 12r_, 
-l + 3r_ +3f+ 




(29) 
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1/3 

1/6 0.25 



- 2/3 




FIG. 12. Section of the intersecting tetrahedron with the 
separating one-leaf hyperboloid. 



Assuming that a convex decomposition along (|28|) is 
possible, we thus arrive at a threefold case distinction: 

• The line segment lies entirely in M\. 

Then it must be tangent to the hyperboloid M*, 
and also an affine direction for R\. The vector t is 
uniquely determined up to a factor by these condi- 
tions. However, that does not mean that the cor- 
responding line segment lies in M\, and, in fact, 
one can show that it never does. Hence this case is 
ruled out. 

• The line segment lies entirely in M 2 . 
This is ruled out analogously. 

• The line segment crosses from Mi into Mi. 
Then t must be affine for both radius functions. 
Again, this determines t to within a factor. But for 
a proper decomposition we must have also that the 
slopes of y/Ri and \[R2 match at u — 0. One can 
show that this never happens inside the tetrahedron 
we discuss, so this case is also ruled out. 



We conclude that no point on M* allows a convex de- 
composition inside V\, and the theorem is proved. 
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APPENDIX: ANALYSIS OF P 

In this appendix we give a characterization of the sep- 
arability classes (T, B\ , and Vy) of P showing that they 
can be deduced from those of P without any computa- 
tion. 

The intimate relation between the two twirls emerged 
already in Lemma ^ where we stated the existence of 
an isomorphism between the two algebras spanning the 
eigenspaces of P and P. This isomorphism establishes 
an affine mapping i between the two eigenspaces that 
we used to compute V\. Due to the inclusion T C 
B\ C V\ it is clear that the same mapping transports 
the sets T and B\ to their counterparts T and B\. 
The mapping t can be computed by fixing the order- 
ing {1, X, V, VXV, XV, VX} for the second algebra and 
concatenating the transformations [| and |2^ getting 



with 



/ d- 1 
I d+1 









d+2 



2d+2 
d-2 



d+1 

d-1 2d-2 

2 2 1 

d+1 d-1 d?-l 

2 2 d 
d+1 d-1 





d+2 
2d+2 
2-d 



V 



o 











2d- 2 
d 

1 

-l d^T 













o \ 







^ I 



With this mapping we can compute directly the P- 
projection of the states A to G: 



(1/2,1/2,0,0,0) 
2 

0, -, — -,0,0 



A : |123) 

C : (|111) - v/3|H2) + V3|121) - 3|122»/4 

'4 + 5d 3d -6 d + 2 



D : 1 122) - 
£:(|112) 

F : (|123) 



> (1,0,0,0,0) 
|121»/V2- - ! (I 
|132))/V2- 



8d - 8 ' 4 - Ad 2 
d-2 1 



0,0 



d- 1' 1-d 
(0,1,0,0,0) 



0,0 
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G : (|112) - |121> 



V3|122))/V5 
'3 2d-A 2 
5'5d-5' 5-5cf 



0.3 



0.35 



0,0 



Applying the transformation to the extremal points 
and inequalities of Theorems [j| [?] and [l2] yields then a 
characterization of T , B\ and V\ . 

We omit here the results of these transformations and 
give the picture corresponding to Fig.[| 



0.50 



0.25 




0.25 0.50 l 

FIG. 13. Sections and projections of T and Si with/onto 
the s+-s_-plane for d = 3. Black: section with T, dark grey: 
projection of T, light grey: section with B\. 



In contrast to what can be seen in Fig. || the projection 
of T onto the s+-s_-plane differs from its section with it 
as one can see in Fig. 14 : 



0.25 




0.25 



0.2 



0.2 



0.3 0.35 
FIG. 14. Zoomed region of Fig. 11. 
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